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Summary 

The theory of elasticity evolved over centuries through the contributions of eminent scientists like Cauchy, Navier. Hooke 
Saint Venant. and others. It was deemed complete when Saint Venant provided the strain formulation in 1860. However unlike 
Cauchy, who addressed equilibrium in the field and on the boundary, the strain formulation was confined only to the field. 

Saint Venant overlooked the compatibility on the boundary. Because of this deficiency, a direct stress formulation could not be 
developed. Stress with traditional methods must be recovered by backcalculation: differentiating either the displacement or t e 
stress function. We have addressed the compatibility on the boundary. Augmentation of these conditions has completed the 
stress formulation in elasticity, opening up a way for a direct determination of stress without the intermediate step ol 
calculating the displacement or the stress function. This Completed Beltrami-Michell Formulation (CBMF) can be specialized 
to derive the traditional methods, but the reverse is not possible. Elasticity solutions must be verified for the compliance of the 
new equation because the boundary compatibility conditions expressed in terms of displacement are not trivially satisfied. This 
paper presents the variational derivation of the stress formulation, illustrates the method, examines attributes and benefits, and 
outlines the future course of research. 


1.0 Introduction 

The theory of elasticity evolved over centuries through the contributions of eminent scientists like Cauchy, Navier. Hooke, 
Saint Venant, and others (see fig. 1 ). Saint Venant. through his strain formulation, presumed to have provided the last set of 
equations in 1860. Decades earlier, Cauchy had formulated the equilibrium in the field and on the boundary of an elastic 
continuum. The equilibrium and compatibility formulations were deemed complete even before the turn of the 20th century. 
This presumption is in error. Saint Venant, credited for the field equations, overlooked the boundary formulation. In other 
words, the theory of elasticity camouflaged a deficiency in the compatibility formulation for well over a century. This weakness 
diverted the development of a direct stress determination method. The traditional methods recovered stress by backcalculations. 
either by differentiating the displacement or the stress function. We have addressed the boundary compatibility condition 
(BCC). Augmenting available elasticity equations with the BCC has led to the completion of the Beltrami-Michell formulation 
(CBMF). This method with stress as the primary unknown can solve displacement and mixed boundary value problems in 
elasticity. A direct stress determination method, like the CBMF, however was envisioned by Beltrami and Michell (ref. 1 ) 
following the strain formulation in 1860. Michell’ s thought is expressed by Love (ref. 2) in the following quotation: 

“It is possible by taking account of these relations [the compatibility conditions ] to obtain a complete system of 
equations which must be satisfied by stress components, and thus the way is open for a direct determination of stress 
without the intermediate steps of forming and solving differential equations to determine the components of 
displacements.” 

Their method, or the classical Beltrami-Michell formulation, was incomplete because it did not include our boundary 
conditions. It thus has limited application. The completed method, or the CBMF, is the versatile elasticity formulation. The 
CBMF can be specialized to obtain Navier s displacement method and Airy's stress function formulation. The CBMF cannot 
be derived from the later two methods. Elasticity solutions that have been obtained by the two traditional methods must be 
verified for the compliance of the new equations because the boundary compatibility conditions expressed in terms of the 
displacement functions are not trivially satisfied. This paper presents the CBMF for three-dimensional elasticity in the 
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subsequent sections: completed Beltrami -Michel I stress formulation, properties of the compatibility conditions, illustrative 
examples, a discussion on attributes and benefits, future course of research, and conclusions. 
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Figure 1 . — The theory of structures had camouflaged a deficiency in the compatiblity formulation since 1860. 


NASA/TP — 200 1-210515 


i 


2.0 Completed Beltrami-Michell Stress Formulation 

We emphasize that the CBMF can be derived only through a variational approach. The new equations, the boundary 
compatibility conditions, are not amenable to a direct derivation. The CBMF is obtained from the stationary' condition of the 
variational functional (ref. 3) of the Integrated Force Method (IFM). The stationary condition yields the field equations and the 
boundary conditions as well as the displacement continuity conditions. The functional n s has three terms— A. B, and W— as 
follows: 


7t v = A + B-W 


(1) 


Term A represents the internal energy expressed in terms of stress (c) and displacement (/*, v, u\) as follows: 


, , du dv dvr [ du dv ^ 

Am l 


dv dw \ 


i du dw 

+, “lfc + aT^ v 


(2a) 


Term B represents the complementary internal energy. It is expressed in strain (£) and the stress function (<p) as 
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(2b) 


The term W represents the potential of the work done due to the prescribed traction on surface P. the prescribed 
displacements on P, and the specified body forces in the volume V of the continuum (see fig. 2) as 


W = 


J (P v H + P v v + P-w)rfS+ | ( P,it + P V V + P- \\ )ds + j" ( B x u + B y v + B. ir) dV 
V s r" V 


(20 


where a„ a v . o ; . x„, T v; . and x. v are the six stress components; £,, e v . £ ; , y„, y r; . and y*arethe six strain components; u. y. and w 
are the three displacement components; cpi, <p 2 > and l P ; are the three stress functions; P ,. P , . and P are the three prescribed 
tractions; and P are reactions where displacements 17 , v . and vv are prescribed; and B,. B, . and K are the three body 

force components. 

The three stress functions are defined as 


9~(p3 5‘(p-> 
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dy 1 9; 2 
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The stationary condition of x s with respect to the displacements and stress functions yields 
as follows; 


a 2 <Pi 

dy 2 

(3) 

all the equations of the CBMF 


NAS A/TP— 200 1-210515 


3 



8tt v = 


Jda, 5 x.xt 3t-v _ L f^T.w „ l5 

{ a.v 3v dz I dv dv 3; 1 


dt- v dx do. 

—^- + —— + —--+6. }■ 6vr + 
dv dv d~ 


3 ' e v d-£- ^'Yvr _ 


d~e v d~£- 3"y v - 

— + ^-iAh-- + 


d 2 e v 3 2 £ v d Yxv 

i-+. — i- 


dv 2 a.v 2 dA-dv 


Stf>3 


dV + 


/[{■ V.Y^.V + «vv t vv +^v;' t iv ~ P, js« +{«v.v T .n +«vy°v +fl v: T y: ~ P v }&’ + 


{ n v.v X 


+ fl vv T v . +« v -<7- - P.}5 h’-|— 



a vz £ \ - — Y>- l+^l flvv®: ^ Yy; |f&Pl - 


< — Y-. + 

dv v ' - 2 ' ac 1 


Idyl " Vvt ' T 7 '' + dv 


dS + 


«v.-e.v — 7 -Y.v; jj&Pi - 

Ov.v«y J-Y.n- j| 5< P3 

J [(« - u )6(o vv ct v + fl vv X xy + a v .X. x ) + (v - v )5(a vv T vv , + a vy a y + a v .x y . )+ 
(vr- i?)5(o V-f T. v +o vv t v . +a vv O- )J dS 


(4) 


2.1 Field Equilibrium Equations 

The three equilibrium equations in the field are obtained as the coefficients of the three incremental displacement 
components (8w, 8v, 5 w) in eq. (4). 


do.. dt n 
— - + — — +• 
d.v dy 

dt XT da v 

— — + -tr^ + 
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ax- r dT v . 
a.v av 


dl- v 
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dz 


-+B- = 0 


(5) 


The formulation of the field equilibrium equations (5) is credited to Cauchy (see fig. 1). 
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2.2 Field Compatibility Conditions 

The three CC in the field are obtained as the coefficients of the three incremental stress (see fig. 1 ) functions (59,. 5<P:- 
5(p,) in eq. (4). The strain CC in the field, referred to as the “strain formulation." is credited to Saint Venant. For an isotropic 
material with Young's modulus E and Poisson's ratio V. the CC expressed in terms of stress are as follows: 



-VO- -VO, 



^2 X 



-vo A 



~VOy - VO- )-2(l +v) 


d-T-v 

dzdx 



-vo v 


-VO- )+— T (Oy -VO- -vo, )-2(I + v) 

<).v _ 


5-T A y 

d.xdv 


= 0 


( 6 ) 


The field CC given by eq. (6) can be further simplified to obtain the canonical form given in reference 4. 

2.3 Boundary Equilibrium or Traction Conditions 

The three equilibrium equations on the boundary are obtained as the coefficients of the variational displacements (8„, 8 V , 
8j in the surface integral terms in eq. (4): 

*v.v°.v + a w^.vv + a vz x .\: ~ P x 
a x>x x xy + a vy o y + a v: t v - = P y 

X V .a vv + Ty-flyy +Oy a V - = Py (7) 

The traction conditions credited to Cauchy are indeterminate because six stresses are expressed in terms of three equations 

2.4 Boundary Compatibility Conditions 

The BCC are obtained as the coefficients of the variational terms (8<pi, S<p 2 , StpO in the surface integrals in eq. (4). For an 
isotropic material, the BCC when expressed in stress become 

j-{o v - (o v -VO- -vo,)-a V y ( 1+v ) T v;} + J“{ a vv ( CT ; ~ VCT .V -V0y)-a v - (l + v)x v; } = 0 

|-{o vt (a. -vo, -VOy )— a v . (l+v)x.,}+|-{a v - (o, -vo, -vo-)-a v , (l + v)x.,} = 0 (8) 

^-{n vv (a, -vo, -vo.)-a vv (l+v)x,y }+^-{a v ., (a, -vo. -vo,)-a vv (l + v)x, v } = 0 

The three BCC given by eq. (8) had been missing since 1860. The traction condition given by equation (6) and the BCC 
ensure stress determinacy on the boundary . 
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2.5 Displacement Continuity Conditions 

Displacement boundary conditions are obtained as the coefficients of the variational reaction terms { 5(fl vv <7 v + a vy T xy + 
(iv-T-x)* 8(a vx r vv + ...), and <S(a vv T- v + ...)} in the surface integral in eq. (4); 


it = u 

v = V (9) 

u' = vT 

The variational functional of IFM yields almost all the equations of elasticity: 

( 1 ) The EE in the field and on the boundary: These are identical to Cauchy’s equations. 

(2) The CC in the field and on the boundary: The field equations are credited to Saint Venant. The BCC are the new 
equations. 

(3 ) The displacement continuity conditions: The variational method provides a difficult, but elegant, derivation of the 
continuity conditions. 

(4) The IFM functional can be specialized to obtain the potential and the complementary energy functionals. 

CBMF is defined through the EE and the CC in the field and on the boundary of an elastic continuum, given by eqs. (5) to 
(8). Stress determination by CBMF does not require the displacement boundary conditions. The displacement functions, u( a*, y, 
:), v(a\ y, :), and u (a, y, c), if required, can be determined by integrating the known stress. The calculation of the constants of 
integration in the displacement functions requires the displacement boundary conditions given by eq. (9). CBMF can be used 
for analysis of stress, displacement, and mixed boundary value problems. A correct solution must satisfy all elasticity equations 
(eqs. (5) to (9)). Solutions obtained without the use of the BCC should be verified for their compliance, especially for a 
continuum with deformable boundaries. In traditional solutions to the elasticity problem, displacement continuity conditions are 
augmented with the slope or derivative of displacement. The ‘'slope" condition somewhat resembles the boundary compatibility 
condition. 


3.0 Properties of the Compatibility Conditions 


Strain or deformation balance is the physical concept behind the controller-type compatibility conditions. In elasticity, 

for example, the strains (e) are controlled,/ (e v , e v , .... y^) = 0, or the deformations <p) are balanced,/ (p b f$ 2 (} n ) = 0, in 

discrete structural systems. The concept of balancing stress and force on a differential block to generate the EE cannot be 
applied to generate the CC. The procedure adopted by Sokolnikoff (ref. 4) — the elimination of displacement from the strain 
displacement relations to derive the field CC — cannot be extended to generate the BCC. Their generation requires the 
utilization of variational calculus. The BCC derived from the stationary condition of the IFM variational functional can, 
however, be verified through an application of the theorems of integral calculus. 


3.1 Verification of the Boundary Compatibility Conditions 

Green's theorem can be used to verify the BCC. This theorem in two-dimensions can be written as 



( 10 ) 


where P and Q are differentiable functions and T represents the boundary of the domain S . 
The BCC in two-dimensional elasticity can be written as 
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This condition is recovered in the line integral term in eq. (10) when the field compatibility condition given by 


3e v 3e v 


_ a 

y 

*y ] 

1 

/ 

ay.n 


a 

/ 

ae.v 

Vi 

ay.nj 

3y 2 ’ a.v 2 

dxdv 

~dx 

A 

a.v 

2 

V 

dy 

> 

a.v 

A 

dy 

J 2 

a.v 

/j 


is substituted in the surface integral. The reader can verify the BCC for three-dimensional elasticity by repeating the proceduie. 

3.2 Nontrivial Property of the Boundary Compatibility Conditions 

The compatibility conditions in the field are automatically satisfied when expressed in terms of continuous displacement 
functions («. v. u). leading to an identity like {/(«. v. n») -flu. v, »•)} = 0. The boundary' compatibility conditions expressed in 
terms of continuous displacements are not satisfied automatically. For three-dimensional elasticity, the three BCC written in 
terms of displacement take the following forms: 
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For three-dimensional elasticity, the compatibility at a boundary interface, defined by an outward normal (n) with the 
direction cosines a, y . a n , and a,. . shown in figure 2. will be satisfied provided the following three residues are matched at the 
left and right of the boundary'. 


% L + 9t,* = 0 *2* + = 0 91/ + 9v/ = 0 (12) 

where 91/ and 91/ refer to the residues for the displacement functions on the right and left of the boundary P-Q in figure 2. 
The compatibility compliance imposes a restriction on the derivatives of the displacement functions. 


NAS A/TP— 2001 -210515 


7 



Figure 2. — A composite elastic continuum. 


3*3 Compliance of the Boundary Compatibility Conditions 

The compliance of the BCC for membrane analysis is illustrated considering an example of a simple stiffness-method- 
based finite element model with a four-node rectangular and a three-node triangular element as shown in figure 3. Along the 
interface connecting nodes (2 and 4) the compatibility compliance is defined by 



1 2 

Figure 3. — Boundary compatibility compliance for 
two-element model. 


where the residue function (9t) for each element defined in terms of displacements (n, v) is 



\ a 2 * if a 2 v a 2 * j 

dxdy 2 fa 2 dxdv I v ' dxdv 2 ^ v - dxdv 
L V /J L /_ 


Consider displacement functions for the rectangular membrane elements (ref. 5) as follows: 


rectangular 


w(-V. y) = C|.y + c 2 xy + cyv + c 4 


(14) 


(15) 
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Likewise, the rec,an * ular v>(.Y, y) displacement function can be defined. The eight constants (cj. c 2 . .... t’s) can be linked to the eight 
nodal displacements of the rectangular element. 

Displacement function for the triangular element can be written as 


lriansular u(.v. y) = d |.v + d 2 y + di 

Likewise the ( ,rianfula 'v(.v, y» displacement can be written. The six constants (</,. d 2 d 6 ) can be linked to the six nodal 

displacements of the triangular element. 


The contribution to the residue function (1H) for each of the two elements can be obtained as 


^a-canjlan >--> = 0.5 { Cl n C b + a u C 2 } 

^triangular! 1-2) _ q (17) 

The boundary compliance (iR rcctal ’ sula, + 9T angular = 0) at the interface of the finite element model shown in figure 3 is not 
satisfied because 9T lcr,a “ = 0.5{a iv c 6 + a vv c 2 } and c 6 = c 2 * 0. In finite element analysis, the traditional assumption that the 
stiffness method satisfies the compatibility conditions a priori needs to be reviewed with respect to compliance oi the boundary 
compatibility conditions. 

3.4 Attributes of Elasticity Operators 

The equilibrium operator and the compatibility operator are related. The operation of the field equilibrium operator [£ e ] on 
the field compatibility operator [Cq] produces a null condition: 


lAiKf ={o) 

The field equilibrium equations, in operator form, can be written as 


dx dv 

0 - I- 

dx dv 

Likewise, the strain compatibility conditions in the field can be represented as 
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(18a) 


(18b) 


(18c) 
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Equation (18a) can be verified by direct substitution. The equilibrium and strain displacement operators also are the 
transpose of each other. 


tel =14 


scir 


(18d) 


Jx 


— 0 


fv 

lY.n 


0 — 

d_ 

dv 


— [Adi* 3 


The strain displacement relation retains the genus of the equilibrium and compatibility concepts. Both the field CC and the 
field EE can be derived from the strain displacement relations by eliminating the displacements and applying virtual work 
principles, respectively. The null property given by equations (18a) and (18b) shows the relationships that exist between the 
equilibrium, compatibility, and strain displacement concepts. The equilibrium and the compatibility concept remain the same 
either in the field or on the boundary. Their equation forms differ to accommodate the field or boundary. 


4.0 Illustrative Example: Composite Circular Plate Under Thermomechanical Loads 

The CBMF is illustrated through the analysis of a composite plate. The moment and displacement of the problem are 
functions of a single radial coordinate. This trivial problem, solved in closed form, has identical solutions via CBMF and 
Navier's displacement method, because for the radial symmetric plate the boundary compatibility condition reduces to the 
slope continuity condition. The problem, however, illustrates the treatment of the boundary compatibility conditon and the 
displacement boundary condition in the stress formulation. 

The clamped composite circular plate is made of an aluminum inner plate (£3,) and a steel outer plate (Q„) as depicted in 
figure 4. The inner plate carries a uniformly distributed load of intensity q , whereas the outer sector is uniformly heated with a 
temperature differential of AT. The radii of the inner a and outer b sectors are 6 and 12 in.; the thicknesses of h, and h () are 0.2 
and 0. 15 in.; the Young's moduli E\ and E 0 are 10.6 and 30.0 million psi; the Poisson's ratios v, and v„ are 0.33 and 0.3; and 
the load q and temperature differential AT are 100 psi and 50 °F, respectively. 

The classical Beltrami-Michell Formulation cannot be used to solve the problem because the plate has clamped, or 
displacement boundary' conditions, at the outer edge and a mixed condition at the interface. CBMF, however, can solve 
the mixed boundary value problem. CBMF considers the two moments (Af r and M<p) as the primary unknowns. The single 
displacement function (w) is obtained from the moments by backcalculations. The CBMF formulation for the problem is 
obtained from the stationary condition of the IFM variational functional. Only the final equations are given. 
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Figure 4.— Composite circular plate subjected to uniform load q and 
temperature differential At. (a) Side view of plate, (b) Planform of plate. 


4.1 Field Equations 

In the field the problem has one equilibrium equation and one compatibility condition. The two equations ensure the 
determinacy of the moment in the field. The field EE is 



dMy 

dr 


+ rq = 0 


The field CC is 

r ^:K- vM r)+(i+v)K-M r )+Ar/-^-|^-j=o 
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where the plate rigidity is K = Eh7[ 1 2( 1 — v : )]. 

Field EE ( 19) and (20) apply to both plate segments (Q„ and Q,)* 

4.2 Boundary Conditions 

The problem has three boundary conditions, consisting of two equilibrium equations and one compatibility condition. At 
the outer clamped boundary the transverse and rotational equilibrium conditions yield two equations: 


R »=-J ~(rM r )-M v ( 21 ) 

dr T 

M, r ' 4 — ' 

where R u and R X{ represent the transverse and moment reactions. The two boundary equilibrium conditions given by 

equations (21) and (22) are used to back calculate the two reactions at the outer boundary but are not explicitly used to calculate 
the moment functions. 

The single boundary compatibility condition at the outer boundary has the following form: 


_ 1 _ 

K 




(23 BC X ) 


This condition, henceforth referred to as BC } will be used to determine the moment functions in the CBMF. 

At the interface, there are three boundary conditions representing the residue balance condition for equations (21) to (23), 
discussed earlier. The three conditions at the interface (r = a) are 


Ml = A/, n 


(2450,) 


d_ 

dr 


{rM\)-M\ = 



(25 «C 3 ) 


-L«-vX) +a ,'5 = -^K-v"«?) + a,"^ 


(2 6BC 4 ) 


where the superscripts I and II refer to the two sectors at the interface. 

The residue in the rotational and transverse equilibrium is accounted for through equations (24) and (25), respectively. 
The compatibility residue balance is represented by equation (26). The temperature effect is introduced through the BCC, or 
equation (26). The CBMF Hill use these three conditions, which will be referred to as BC\, BC2, BC 3, and BC4. 

The finite condition of the moment functions (M r and M 9 ) referred to as (BC 5 and BC 6 ) will also be used in the CBMF 
analysis. The two field equations and the six boundary conditions are sufficient for the calculation of the two moment 
functions. In CBMF, displacement is backcalculated from moment functions. This process requires displacement continuity 
conditions (vr = 0 at outer boundary and at the interface: vv = vv°) that are obtained from the stationary condition of the 
variational functional. 
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4.3. Solution 

sThe field equations (eqs. ( 19) and (20)) are solved to obtain the following general solutions. For the inner plate. 


M'r (i ) = --T + ^ c \ 0 + v, )logr + ^C, (1 -Vj)+^A ~(3 + V/ )^- 
<('-) = -4 + T C ' O + v^logr-ic, (1 _ V| )+-^£ > l - ~0 + 3v; )q> - 
Likewise, the solutions for the outer plate with no distributed load are 


M r ( r ) = -~T+\ C 2 (l + V„)log /' + - J Cl (l-V f) ) + ^£>2 


(27a) 


(27b) 


(28a) 


(r) = -- C 2 (l + v 0 )logr-jC 2 (l-v „)+^£> 2 1 

r~ 4 

The six integration constants are evaluated for the six boundary conditions for a temperature AT of 50 °F, a load q of 
100 psi, and the material properties specified earlier. The moment solutions follow. 

For the inner plate (0 < r < b ) 


M' r (r) = 844.05-20.81r 2 (29a) 

A/' (r) = 844.05 - 20.8 Ir 2 (29b) 


For the outer plate (a <r<b). 


M° (r)= 2046.63 - (5203.06 / r 2 )- 1 1701og r (30a) 

Af " (r)= 2676.63 + (5203.06 /r 2 )-11701ogr (30b) 

The moment solutions for the composite plate with the mixed boundary condition are obtained using CBMF without any 
reference to displacements in the field and on the boundary. Displacements, if required, can be backcalculated front the 
moments. For this problem, the displacement conditions at the outer boundary (vv = 0 at r = b) and the continuity conditions at 

the interface (w = w°) are sufficient to calculate the displacement. 

For the inner plate, the transverse displacement (0 < r < a) is 


vv' (r) = 3. 1 209 - 0.0356r 2 + 0. 1 757 x 1(T 3 r 4 (31a) 

and for the outer plate (a <r< b). it is 

w°(r) = 5.3614 -0.1344r 2 -0.72961ogr + 0.04417r 2 logr (31b) 

The CBMF solution procedure illustrated here can be extended to analyze displacement and mixed boundary value 
problems in elasticity as well as for plate and shell problems. 
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5.0 Discussions 


The discussion is given under the role of the CC, benefits front the use of the CC. and why the CC weren't formulated 
earlier. 

5.1 Role of the Compatibility Conditions 

The role of the CC can be perceived from an examination of the stress-strain law, which is universal to all analysis 
formulations. Hooke's law (see fig. 1 ). relates stress {a} and strain {£} through a material matrix [9??] as ( {a} = [OT]{e}). The 
stress must satisfy the state of equilibrium, whereas the strain must satisfy the compatibility. Hooke's law, the equilibrium, and 
the compatibility are sufficient for determining the stress in an elastic continuum. Conceptually, the calculation of stress can be 
represented as 


Equilibrium Equations 
Compatibility Conditions 


{Stress} 


Mechanical Load 
Initial Deformation 


(32) 



Figure 5. — Equilibrium equations and compatibility 
conditions in elasticity. 


Displacement in explicit terms is not required for the determination of stress. The equilibrium and compatibility concepts 
of elasticity are depicted in the two halves of the pie diagram in figure 5. The immaturity in the compatibility condition at the 
boundary’ is represented by the shaded quarter. In structural mechanics, the compatibility formulated through the concept of 
redundant forces, "‘cutting" and “closing" the gap, is not even parallel to the strain formulation of Saint Venant. In other words, 
the pie diagram also applies to structural mechanics. The theory of elasticity has been developed utilizing the information 
contained in the three-quarters of the pie diagram. The utilization of the additional quarter should improve the theory making it 
more robust. 

Stress can be determined bypassing the compatibility conditions, but omitting the compatibility can lead to erroneous 
stress. Even a century ago, this deficiency was observed by Todhunter (ref. 6) while he was scrutinizing astronomer Royal 
Airy's (1801-1892) attempt to analyze the stress in cantilever and simply supported beams. 

“Important Addition and Correction. The solution of the problems suggested in the last two Articles were given — 
as has already been stated — on the authority of a paper by the late Astronomer Royal, published in a report of the 
British Association. I now observe, however — when the printing of the articles and engraving of the Figures is already 
completed — that they cannot be accepted as true solutions, inasmuch as they do not satisfy the general equations (164) 
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of § 303 [note that the equations in question are the CC ]. It is perhaps as well that they should be preserved as a 
warning to the students against the insidious and comparatively rare error of choosing a solution which satisfies 
completely all the boundary conditions, without satisfying the fundamental condition of strain [note that the condition 
in question is the compatibility condition], and which is therefore of course not a solution at all.” 

A solution to a problem can be obtained utilizing a subset of elasticity equations. Such a solution was obtained by Airy. Its 
inaccuracy was pointed out by Todhunter. Strictly speaking, a valid elasticity solution must satisfy the boundary compatibility 
conditions. It is quite possible that many solutions obtained earlier may satisfy the new conditions without explicit imposition. 
The “process of evolution” might have eliminated inaccuracy in the traditional solutions because such results have been in 
existence for over a century. Despite this conciliatory concession, it is but prudent to verify the traditional solution for the 
compliance of the new boundary conditions. 

5.2 Benefit in Finite Element Analysis 

A finite element method that parallels the completed Beltrami-Michell formulation has been developed. This is referred to 
as the Integrated Force Method, or IFM. Force parameters are its primal variables (see table I). Displacements are recovered 
from forces. The IFM equations for a finite element model with n force and m displacement unknowns are obtained by coupling 
the m equilibrium equations ([B]{FJ = (F()to tbcr=n-m compatibility conditions ([C]|G]{ F] = {§/?)): 


[B] 

[C][G] 



or [S]{F} = {/>*} 


From the internal forces {F} the displacements {X} are backcalculated as 


M = [J]{[g]{f}+{p 0 }} 


(33) 


(34) 


where [J] = m rows of { [S] ! } 1 . 


TABLE I.— METHODS OF STRUCTURAL MECHANICS AND ASSOCIATED VARIATIONAL FUNCTIONALS 


Method 

number 

Method 

Primary 

variables 

Variational 

functional 


Elasticity 

Structures 

Elasticity 

Structures 


1 

Completed Beltrami-Michell 
Formulation (CBMF) 

Integrated Force 
Method (IFM) 

Stresses 

Forces 

IFM variational 
functional 

2 

Airy's formulation 

Redundant force method 

Stress function 

Redundants 

Complementary 

energy 

3 

Navier's formulation (NF) 

Stiffness method (DM) 

Displacements 

Deflections 

Potential energy 

4 

Hybrid method (HF) 

Reissner’s method (RM } 

Stresses and 
displacements 

Forces and deflections 

Reissner’s 

functional 

5 

Total formulation (TF) 

Washizu's method (WM) 

Stresses, strains, and 
displacements 

Forces, deformations, 
and deflections 

Washizu’s 

functional 


In equations (33) and (34), [B] is the m x n rectangular equilibrium matrix, [G] is the n x n flexibility matrix, [C] is the 
rxn compatibility matrix, {bR} = -[C]((3 } is the r-component effective initial deformation vector, ( P ) is the initial 
deformation vector of dimension /j, [S] is the IFM governing matrix, and [J] is the m x n deformation coefficient matrix. 

A research-level finite element code, referred to as IFM/Analyzers, has been developed. Results for a few r examples 
obtained using the IFM/Analyzers code and a commercial stiffness method code are depicted in figure 6. The three IFM 
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Figure 6. — Solution to a set of examples. 






















elements used are QUAD0405, which is a four-node membrane element with five force unknowns; HEX2090. which is a 
20-node brick element with 90 force unknowns; and PLB0409. which is a four-node plate bending element with nine unknown 
forces. These are very simple elements that use standard interpolation functions. The elemental generation used numerical 
integration but did not use reduced integration or bubble function techniques. 

Displacement and moment solutions for a flat plate under a concentrated load is given as figure 6(a). A displacement 
solution to a cantilever beam using membrane elements is shown in figure 6(b). A stress solution to a cantilever beam modeled 
with brick elements is shown in figure 6(c). The frequency analysis of a turboprop blade using a brick element is shown in 
figure 6(d) reported under problem 5. The test results and stiffness method predictions are shown in figure 6(e) for a beryllium 
cylinder. For the examples the integrated force method outperformed the stiffness method, overshadowing the simplicity at its 
element level. The developer of the commercial code disputed the IFM results (ref. 7). To resolve the dispute, the examples 
were solved again and a detailed report with several tables was prepared. No numerical error could be detected in any of the 
example problems. Only one typographical error was found— in the element subtitle. The correct word is “Six/HX20_90.” not 
"Three/H X08_90' ' as published in the journal (ref. 7). The information submitted to the journal and published in the NASA 
Technical Memorandum (ref. 8) was correct. The typographical error in the publication is inconsequential to the conclusion. 

5.3 Why Weren’t the Compatibility Conditions Developed Earlier? 

Three reasons for the tardy development of the compatibility might be complexity, complacency, and emphasis on 
solutions. Compatibility conditions in relative terms can be considered to be more complex than other elasticity relations 
because their formulation requires the variational concept of calculus. If compatibility were simpler, then its derivation would 
have matured earlier. However, the converse would not be a certainty. In figure I, we have depicted a partial list of pioneer 
scientists who contributed to the development of the subject. Complexity cannot be considered to be an insurmountable 
obstacle with Cauchy. Saint-Venant, Bernoulli, or Maxwell, to mention just a few. Complexity is not the likely reason. 

An argument in favor of complacency cannot be justified. Since eminent scientists developed the basic principles of this 
science, they could not have overlooked the fundamental compatibility concept. In the theory' of structure, consider the 
displacement and classical force methods. Displacement {X} is the unknown in the displacement method [K]{X) = {P). 
However, the redundant { /?}, not the internal force {F}. is the unknown in the classical force method. [H ] { = { P). 

Since displacement is the unknow n of the displacement method, force should have been the unknow n of the force method. 

In elasticity, Navier's displacement method can solve all three types of boundary value problems. The classical 
Beltrami-Michelfs stress formulation cannot be used to solve the displacement or mixed boundary value problems. Even 
novice researchers would surely have noticed and questioned the lack of uniformity between the force and displacement 

methods. . . 

Even an approximate solution can fulfill the need of an industry, which may have little interest in the particular method 

used. The building and bridge industries required analysis of indeterminate trusses, continuous beams, frames, and other 
structures. Redundant analysis solved such problems manually. For skeletal frames, the method of moment distributions and 
Kani's method of stiffness balance were very popular prior to computer automation. Elasticity results, obtained through 
superposition techniques have served the industry well. The results required by civil, mechanical, and aerospace industries were 
given primary importance, but the analysis methods used were considered but secondary. The point is that the solutions of 
engineering problems not only became central to the work but occupied most of the available time of competent researchers, 
leaving little or no time for them to ponder or address the deficiencies or completion of the theory of compatibility. The 
engineering problem-solving aspect of structural mechanics was considered to be most important. Such jobs were considered to 
be glamorous and paid high dividends. But the basic science of structural mechanics, including the completion of the theory of 
compatibility, appears to have been neglected by industry', research institutions, and academia alike. 

Glamour and the many other dividends associated with solving difficult industrial problems are. in our opinion, the 
primary reasons behind the slow progress of the theory of compatibility. Complexity and complacency may be considered but 
secondary reasons for neglecting such conditions. 
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6.0 Future Research 


The method of force or the completed Beltrami MichelFs formulation and the method of displacement or Navier's method 
are the two primary formulation of elasticity. Other methods are listed in table 1. The boundary compatibility allows facile 
movement between the five methods given in table 1: (1 ) CBMF, (2) Airy's stress function method. (3) Navier's displacement 
method. (4) Reissner's method, and (5) Washizu's formulation. It is prudent to obtain solutions to elasticity problems by both 
force and displacement methods and to eliminate error by comparison. Research thus far emphasized Navier's displacing 
method, which might have entered the plateau of diminishing marginal return. The force method has opened up, and 
researchers should exploit its potential. We propose initial research under two subtopics: completeness of the theory and 
verification of available solutions. 

6.1 Completeness of the Theory 

We must complete the theory of linear and nonlinear elasticity. We have addressed the theory of linear elasticity. The IFM 
variational functional has been formulated (ref. 3 and 9). and the BCC have been generated for two- and three-dimensional 
elasticity problems in Cartesian coordinates. The BCC were also derived for two-dimensional problems in polar coordinates, 
which has yet to be published. The BCC have been generated for a rectangular plate flexure problem (ref. 10), a circular plate, 
and a radially symmetrical cylindrical shell (ref. 9). Solutions to plate and cylindrical shell problems have been obtained for 
mechanical load and temperature distribution. CBMF has been extended to vibration and buckling problems utilizing the 
concept of “stress mode shapes*' (ref. 11). In addition, the noncompliance of the CC at the boundary has been shown for a 
simple elasticity problem (ref. 1 1 ). CBMF needs to be developed in curvilinear coordinates for three-dimensional elasticity. 

The procedure should be extended to plates and shells. Since there are many different geometrical shell configurations (like 
cylinder, spherical, conical, and paraboloidal shapes), the feasibility of training the IFM variational functional through a 
computer program using symbolic language should be explored. The linear CBMF analysis should be extended next to finite 
deformation elasticity theory in Lagrangian and Eulerian coordinates. The steps followed for linear elasticity can be adopted 
for nonlinear problems: an extension of the variational functional including stress functions, variational operations, recovering 
equations, and their interpretation and attributes. Prior to three-dimensional nonlinear analysis, the basic steps can be exercised 
for a two-dimensional model. 

6.2 Verification of Available Solutions 

Solutions are available to many elasticity problems.. We must verify these solutions for the compliance of the BCC. This 
exercise might be trivial for some problems, but new solutions might be required for others. This exercise cannot be avoided 
because classical solutions are still in use to verify numerical solutions to engineering problems. For three-dimensional 
elasticity, the compatibility conditions on a boundary surface will be satisfied provided the following curvature terms vanish: 


d 2 u 

d 2 u 


_ d 2 u 

dy 2 

~dz 2 ~ 

dzd.X 

d.xdy 

d 2 v 

3 2 v 

d 2 v 

_ a 2 v 

dy 2 

~ dz 2 ~ 

dydz 

dxdy 

d 2 w 

_ d 2 w _ 

d 2 w 

_ 3 2 vv 

d.x 2 

dv“ 

dydz. 

" dzdx 


The boundary compatibility conditions can be factorized to obtain: 
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( 36 ) 


Boundary compatibility compliance requires equality ol the following rotation teims. 


5\v _ dv d\v _ dit 

dv dz dx dz 

The displacement w with respect to y (or in the ^ -displacement direction), must be equal to the derivative of the 
displacement v with the respect to c (or vv-displacement direction), and so forth. 


dv _ du 
dx dv 


(37) 


7.0 Conclusions 

A deficiency in the theory of elasticity has been eliminated through the formulation of the compatibility conditions on 
the boundary of an elastic continuum. These equations have completed the stress formulation in elasticity a task that was 
originally attempted by Beltrami and Michell at the end of the 19th century. This primal method of stress can handle all 
three classes of boundary value problems for static, thermal, and dynamic loads. Elasticity solutions generated by Navier s 
displacement method and Airy's stress function formulation should be verified for the compliance of the novel compatibility 
conditions because these equations expressed in displacements do not become trivial functions like the field compatibility 
conditions. Structural analyses and design can benefit from the use of the stress formulation. The compability conditions remain 
to be formulated in nonlinear elasticity. 
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